Ferromagnetism of magnetic semiconductors — Zhang-Rice limit 
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It is suggested that p-d hybridization contributes significantly to the hole binding energy Eb of 
Mn acceptors in III-V compounds, leading in an extreme case to the formation of Zhang-Rice-like 
small magnetic polarons. The model explains both strong increase of Eb and evolution of Mn spin- 
resonance spectrum with the magnitude of valence-band offsets. For such a structure of Mn impurity 
in III-V materials, possible models accounting for the recently determined Curie temperature of 
about 940 K in a compensated Gao.g1Mno.09N are discussed. 
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Recent theoretical suggestions Jj], || that ferromag- 
netic ordering may persist above room temperature in 
(Ga,Mn)N and related materials have triggered consid- 
erable fabrication efforts ||, [|, || , which have culminated 
in the determination of Curie temperature Tq m 940 K 
for a film of a compensated Gao.91Mno.09N frown by 
molecular-beam epitaxy (MBE) || . In order to put this 
result in an appropriate perspective we note that the 
magnetic ion density is about twenty times smaller in 
Gao.91Mno.09N comparing to Co, which has the highest 
known value of Tq = 1400 K. This adds a new dimension 
to capabilities of GaN-based structures, whose impor- 
tance in photonics and high-power electronics has already 
been proven. While further experimental studies of this 
new ferromagnetic system will certainly lead to develop- 
ments that are unforeseen today, it is already important 
to model the physical environment, in which the ferro- 
magnetism is observed as well as to indicate parameters 
that are controllable experimentally, and which account 
for the magnitudes of Tq and remanent magnetization. 

In this work, we demonstrate that the high Tq value 
quoted above is consistent with the expectations of the 
mean-field Zener model of the carrier-mediated ferro- 
magnetism, which was elaborated earlier by us and 
co-workers jl], [?]]. We then examine various assump- 
tions and approximations of the model exploring recent 
developments in the field of magnetic semiconductors 
@, |, |[ |l|, @. In particular, by a detail tracing 
of properties of Mn impurity across the ensemble of III- 
V semiconductors we suggest that Zhang-Rice (ZR) small 
magnetic polarons are formed in GaN:Mn. We then pro- 
pose that derealization of these polarons may drive the 
ferromagnetic transition in p-(Ga,Mn)N. We also note 
that if the d 5 /d 6 level resides in the gap and the con- 
centration of donors is greater than that of Mn, double 
exchange involving the upper Hubbard d-band may ac- 
count for ferromagnetism of n-(Ga,Mn)N. 

Figure 1 presents Tq of wurzite Ga^jMn^N calcu- 
lated according to the model jl], [?J and for material pa- 
rameters discussed in detail previously 0. The data are 
shown for magnetization perpendicular to the c axis of 



the wurzite structure, which is found to be the easy plane 
in (Ga,Mn)N for the displayed values of the Mn and hole 
concentrations. As seen, Tq of 940 K is expected for the 
hole concentrations about four times smaller than the 
Mn concentration x in Gao.91Mno.09N. Such a value of 
Tq is consistent with that deduced from magnetization 
measurements at 0.1 T between 300 and 750 K [||. 

We begin the discussion of various features of the- 
ory 01 by recalling that according to the mean-field 
Zener model Q, the properties of the system consist- 
ing of localized spins and itinerant carriers can be ob- 
tained by minimizing the total free energy with respect 
to the spin magnetization M. It is assumed jl], [?| that 
the Mn ions substitute the cations, and are in the 3d 5 
5 = 5/2 configuration, whereas the holes reside in the 
valence band, which is described by the disorder- free 
k-p Kohn-Luttinger theory for tetrahedrally coordinated 
semiconductors. The p-d exchange between the two sub- 
systems, the Kondo interaction, is incorporated to the 
k-p scheme within the molecular-field and virtual crystal 
approximations, whereas the hole-hole exchange interac- 
tion is treated in terms of Landau's Fermi liquid theory. 
The hole free energy, not ground state energy, is adopted 
for the evaluation of Tq p], Q|, an ingredient significant 
quantitatively in the high temperature regime, such as 
that covered by the data in Fig. 1. 

Figure 2 shows the energetic position of the Mn im- 
purity level in III-V compounds, as evaluated by various 
authors from measurements of optical spectra and ac- 
tivation energy of conductivity. A priori, the Mn atom, 
when substituting a trivalent metal, may assume either of 
two configurations: (i) d 4 or (ii) d 5 plus a weakly bound 
hole (d 5 +h = A ). Accordingly, the experimentally de- 
termined energies correspond to either d 4 /d 5 or A°/A~ 
levels. It appears to be a general consensus that the 
Mn acts as an effective mass acceptor (d 5 +h) in the case 
of antimonides and arsenides. Such a view is supported 
by a relatively small Mn concentrations corresponding 
to the insulator-to-metal transition, which according to 
the Mott criterion n 1 / 3 as = 0.26, points to a relatively 
large extension of the effective Bohr radius as- More- 
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FIG. 1: Curie temperature of wurzite Gai-^Mn^N according 
to the mean-field Zener model H, [?J calculated as a function 
of the hole concentration per Mn ion, xh = p/xN a . 



FIG. 2: Experimental energies of Mn levels in the gap of III-V 
compounds according to Refs. 0, p], [u| El ^ j5| E|, 
in respect to valence-band edges, whose relative positions are 
taken from Ref. 12 ill . 



over, the ESR studies of GaAs:Mn reveal, in addition 
to the well known spectrum of Mn d 5 with the Lande 
factor gg = 2.0, two additional lines corresponding to 
g a = 2.8 and gt, — 6 (2^, which can be described 

quantitatively within the k ■ p scheme for the occupied 
acceptor [^2], ^3|. Here, the presence of a strong antifer- 
romagnetic p-d exchange interaction between the bound 
hole and the Mn d-electrons has to be assumed, so that 
the total momentum of the complex is J = 1. In agree- 
ment with the model, the additional ESR lines, in con- 
trast to the g = 2.0 resonance, are visible only in a nar- 
row range of the Mn concentration (24) , which should be 
greater than the concentration of compensating donors, 
and smaller than that at which acceptor wave functions 
start to overlap and to merge with the valence band. The 
antiferromagnetic coupling within the d 5 +h complex is 
seen in a number of experiments, and has been employed 
to evaluate the p-d exchange integral (UNq w — 1 eV |2f| 
in GaAs:Mn, the value in agreement with that deter- 
mined from interband magnetoabsorption in (Ga,Mn)As 

Importantly, the above scenario is corroborated by re- 
sults of photoemission [ pTj , and x-ray magnetic cir- 
cular dichroism (MCD) studies |^9|, |3(J in metallic or 
nearly metallic Gai_j;Mn x As. The latter point to the d 5 
Mn configuration. The former are not only consistent 
with such a configuration but also lead to the value of 
j3N similar to that quoted above, j3N w -1.2 eV [|7|. 
Furthermore, the photoemission reveals the presence of 
two features in the density of states brought about by 
the Mn constituent: the original Mn 3d 5 states located 
around 4.5 eV below the Fermi energy Ep , and new states 
merging with the valence band in the vicinity of Ep [ p8[ . 
These new states correspond to acceptors, as discussed 
above. They are derived from the valence band by the 



Coulomb field as well as by a local Mn potential that 
leads to a chemical shift in the standard impurity lan- 
guage, or to a valence band offset in the alloy nomencla- 
ture. The fact that transition metal impurities may pro- 
duce both resonant and Coulomb states was also noted 
in the case of CdSe:Sc 

In contrast to antimonides and arsenides, the situation 
is much more intricate in the case of phosphides and ni- 
trides. Here, ESR measurements reveal the presence of 
a line with g = 2.0 only (J, fTi} j32| , |33| , |3^| , which is thus 
assigned to d 5 (A - ) centers ||! [F| p2j [33] . Moreover, ac- 
cording to a detail study carried out for a compensated 
n-type GaP:Mn Q|, the ESR amplitude diminishes un- 
der illumination and, simultaneously, new lines appear, 
a part of which exhibit anisotropy consistent with the d 4 
configuration. This, together with the apparent lack of 
evidence for d 5 +h states, even in p-type materials, ap- 
pears to imply that Mn in the ground state possesses the 
d 4 , not A , electron configuration |Q. This would mean 
that the Mn energy in Fig. 2 for GaP |TJ] and, therefore, 
for GaN [O, [l^] (where the valence band is lower than in 
GaP) corresponds to the d 4 /d 5 , not A°/A~ level. Such 
a view seems to be supported by the ab initio computa- 
tion within the local spin density approximation (LSDA) , 
which points to the presence of the d-states in the gap 
of (Ga,Mn)N |2j. In this situation, the spin-spin inter- 
action would be driven by a double exchange mechanism 
involving hopping of d-electrons |3^] , as in the case of 
colossal magnetoresistance manganites, making the mod- 
eling leading to the results presented in Fig. 1 invalid. 

However, there is a chain of arguments which calls the 
above interpretation into question. First, guided by pho- 
toemission results for II- VI compounds |36| one expects 
that the energy of the d 4 /d 5 level will vary little between 
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arsenides and nitrides. This implies that this level should 
reside in the valence band of GaN despite the 1.8 eV va- 
lence band offset between GaN and GaAs. The resulting 
contradiction with the LSDA findings can be removed 
by noting that in the case of strongly correlated 3d elec- 
trons, a semi-empirical LSDA+U approach is necessary 
to reconcile the computed and photocmission positions 
of states derived from the Mn 3d shell in (Ga,Mn)As 



Second, it is known that the p-d hybridization, in ad- 
dition to producing the exchange interaction, can con- 
tribute to the hole binding energy Eb [[38], |3^]. We take 
the hole wave function as a coherent superposition of 
p-stated of anions adjacent to Mn (38). Assuming the p- 
orbitals to be directed towards the ion Mn we find that 
the T2 state has 30% lower energy than that correspond- 
ing to the mutually parallel p-orbitals. This shows that 
not only Kohn-Luttinger amplitudes from the T point of 
the Brillouin zone are involved. In order to evaluated 
Ef,, a square- well spherical potential U(r) = UQ(b — r) 
is assumed Q , whose depth U is determined by the p-d 
hybridization taking into account the above mentioned 
arrangement of the p-orbitals, 
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Here, the values of (3Nq are taken from Ref. @; A e // is 
the distance of d 4 /d 5 level to the top of the valence band, 
which is evaluated to be 2.7 eV in (Ga,Mn)As |27f] , and is 
assumed to be reduced in other compounds by the corre- 
sponding valence band offsets (Fig. 2), and Ueff — 7 eV is 
the correlation energy of the 3d electrons jp7L |3(| . Finally, 
b/a is the well radius in the units of the lattice constant, 
and should lie in-between the cation-anion and cation- 
cation distance, y/3a /4: < b < a /v2- It turns out that 
in the case of GaN:Mn the hole is bound by Mn, even in 
the absence of the Coulomb potential, Ef, = 1.0 eV for 
b = 0.46a o . This demonstrates rather convincingly that 
a large part of Ef, originates indeed from the p-d interac- 
tion, indicating that the Zhang-Rice (ZR) limit p8|, |3t| 
is reached in these systems. The formation of ZR small 
magnetic polarons invoked here has a number of impor- 
tant consequences which are now discussed. 

We begin by arguing that when the ZR polarons are 
formed, the d 5 and d 5 +h states acquire virtually the 
same Lande factor g — 2.0. Quite generally, the role of 
the intra-atomic spin-orbit splitting A Q diminishes with 
the carrier kinetic energy, and thus with Eb. An up- 
per limit of orbital-momentum corrections to the Mn 
Lande factor, A g — g — gs can thus be evaluated by 
considering non-magnetic acceptors, for which Eb is typ- 
ically smaller than that of Mn. In the wurzite GaN, the 
orbital momentum is quenched by the uniaxial crystal 
field, and the Lande factors of GaN:Mg are close to two, 
g[ = 2.07 and g£ = 1.99 go). Hence, for the d 5 +h state 
in wz-GaN:Mn, where the total spin S p = 2, we obtain 
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A<7 = (gs — Sh)/6 ~ —0.01. In cubic materials, in turn, 
the total momentum J p = 1 and Ag = 3/4 for deep d 5 +h 
, in agreement with the ESR findings for 
p3, 24 . We note, however, that the Jahn- 
Teller effect should be particularly strong in the case of 
ZR polarons. The corresponding lowering of local sym- 
metry will lead to the quenching of orbital momentum, 
especially in compounds with a relatively small value of 
A . Hence, also in the cubic nitrides and phosphites, 
S p = 2 and thus Ag w 0. Interestingly, by exploiting the 
hole interaction with host nuclear spins, it becomes possi- 
ble to tell d 5 and d 5 +h states despite their similar Lande 
factors j34|. In particular, the presence of the hole ex- 
plains a surprising broadening of the Mn hipcrfine sextet 
already at relatively small concentrations of Mn impuri- 
ties, observed in InP:Mn ||, || GaP:Mn f|, and 
GaN:Mn [g as well the detection of P and In ENDOR 
signals in InP:Mn Q. 

Another important consequence of the the ZR polaron 
formation is the shift of the Mott critical concentration 
towards rather high values. According to the known re- 
lation between Eb and as , the critical hole concentra- 
tion is p c = 4 x 10 19 cm~ 3 in (Ga,Mn)As and at least an 
order of magnitude greater in (Ga,Mn)N, if no shallower 
acceptors are present. According to the two-fluid model 
0, corroborated by experimental results |^9|, |4^, pi| , 
only a part of spins orders ferromagnetically in the in- 
sulator phase. The Mn spins are ordered in space re- 
gions, in which weakly localized holes reside |J. This 
can be the case of the ferromagnetic Gao.91Mno.09A" sam- 
ple [[|, whose spontaneous magnetic moment constitutes 
only about 20% of the value expected for x — 0.09 and 
S = 5/2. Remarkably, the disorder in Mn positions seems 
to shift Tq to higher values in such a case presum- 
ably even above those display in Fig. 1. On the other 
hand, corrections to molecular- and mean-field approxi- 
mations tend to lower Tq values [|[ [l0|. While our work 
provides a background for detail studies of these effects, 
we expect that the corresponding corrections will not be 
qualitatively significant. Indeed, the short-range poten- 
tial of Eq. 1 is still relatively weak: our corresponds 
to J pd = 0.054 eVnm 3 and J = 0.5 eV of Refs. [g, 0, 
respectively. 

It might appear from Fig. 1 that a further enhancement 
of Tq would be possible by increasing the hole concen- 
tration p. Actually, however, if p — > xNq (i.e., xt — > 1), 
the antiferromagnetic portion of the carrier-mediated 
(RKKY) interaction will reduce Tq |, g @ and 
eventually drive the system towards a spin-glass phase 
Q. Thus, it seems that the compensation by donors, 
despite introducing an additional disorder, constitutes 
the important ingredient of ferromagnetism in Mn-based 
III-V compounds p5| . 

Finally, we note that if indeed, as we advocate here, the 
d 4 /d 5 state resides in the valence band, the d 5 /d 6 level 
may lie below the bottom of the conduction band. If this 
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is the case, and the concentration of donors is greater 
than that of Mn, double exchange involving the d 6 elec- 
trons may account for ferromagnetism of n-(Ga,Mn)N. 
We note in this context that because of the anomalous 
(spin) Hall effect Q, a possible coexistence of p-type 
and n-type regions in compensated semiconductors fl47| , 
and the presence of strong electric fields due to lack of 
inversion symmetry, the evaluation of the carrier type 
and concentration from transport measurements is by no 
means straightforward in thin films of wz-(Ga,Mn)N. 

In summary, according to the discussion presented in 
this paper, while the Mn substitutional impurities act as 
shallow acceptors characterized by the Bohr radius spac- 
ing over tens of A in antimonides, in the case of lighter 
anions the hole becomes tightly localized. In the extreme 
case, such as GaN:Mn, the Zhang-Rice-like magnetic po- 
larons are formed. At sufficiently high Mn concentra- 
tions, the holes tend to delocalize and can then mediate 
long-range correlation between the Mn ions. In an over- 
compensated material this correlation may be transmit- 
ted by d 6 electrons, if the corresponding upper Hubbard 
band resides below the bottom of the conduction band. 
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